Relevance of complex branch points for partial wave analysis 
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A central issue in hadron spectroscopy is to deduce — and interpret — resonance parameters, 
namely pole positions and residues, from experimental data, for those are the quantities to be 
compared to lattice QCD or model calculations. However, not every structure in the observables 
derives from a resonance pole: the origin might as well be branch points, either located on the real 
axis (when a new channel comprised of stable particles opens) or in the complex plane (when at 
least one of the intermediate particles is unstable). In this paper we demonstrate first the existence 
of such branch points in the complex plane and then show on the example of the 7riV Pn partial 
wave that it is not possible to distinguish the structures induced by the latter from a true pole signal 
based on elastic data alone. 

PACS numbers: 14.20.Gk, 13.75.Gx, 11.80.Gw, 24.10.Eq, 



I. INTRODUCTION 

The second and third resonance region of baryonic ex- 
cited states is currently under intense experimental inves- 
tigation at various laboratories such as ELSA, MAMI, 
or JLab [IHIj]. Many resonances overlap at these ener- 
gies, and usually partial wave analyses in different frame- 
works, such as iC-matrix approaches or dynamical cou- 
pled channel models [oT-f^ are necessary to disentangle 
the resonance content. Furthermore, many resonances 
may couple only weakly to the irN channel, and the in- 
vestigation of different initial and final states in hadronic 
reactions is mandatory [l8| . Also, at these energies multi- 
pion intermediate and final states are becoming increas- 
ingly important and should be included in the analysis of 
the iS-matrix. For the corresponding T-matrix, channels 
with stable particles like rjN induce a branch point at the 
threshold energy (y/s = rrir, + Mn), that may be visible 
as cusps in the amplitude [23|, [24| • 

For effective multi-pion channels with one unstable and 
one stable particle, such as pN , the analytic structure is 
more complicated. In comparison to the branch points 
on the real s axis and the first and second sheet poles, 
the third type of allowed singularities are branch points 
within the complex energy plane. They emerge when 
amongst groups of particles of an at least three-body de- 
cay there exists a strong correlations between two parti- 
cles. For example, a significant fraction of ir + ir~X inter- 
mediate and final states typically goes through the p me- 
son. The resulting line shapes are discussed in Ref. [25| . 
Branch points in the complex plane also emerge in the re- 
cently developed complex-mass scheme for baryonic res- 
onances (2f|. 

Known theoretically for a long time [13, H^|, these 
branch points are present in several modern approaches, 
such as the GWU/SAID analysis the Jiilich [3- 

[l8| and EBAC [20|, approaches, or the Bonn- 
Gatchina [l3[ analysis. It is the goal of this study to 



demonstrate the model-independent character of those 
complex branch points. To do so, we employ general 
properties of the 5-matrix only. In a particular example 
it is then shown that the branch points are of relevance 
in partial wave analyses: if the theoretical partial wave 
does not include them, their absence can easily be sim- 
ulated by resonance poles. This, of course, distorts the 
extracted baryon spectrum. Branch points in the com- 
plex plane are thus important for the reliable extraction 
of resonance parameters. 

The paper is organized as follows: in Sec. |TT] the ex- 
istence of branch points in the complex plane is derived 
from three-body phase space, in Sec. Ill Al the properties 
of the branch points are determined, and in Sec. IIIII it 
is shown that these branch points are relevant in the ex- 
traction of the resonance content of partial waves. 



II. ANALYTIC STRUCTURE OF THE 
5-MATRIX AND COMPLEX BRANCH POINTS 

Every channel opening introduces a new branch point 
and with it a new sheet to the S'-matrix, located at s = 
E TO i) 2 i with m i the masses of the stable particles in 
that channel. The first sheet is always the physical one, 
i.e. where the physical amplitude is situated. The only 
singularities allowed on the first sheet are poles on the 
real s axis below the lowest threshold (=bound states) 
or branch points on the real axis. On other sheets, poles 
and branch points can be located anywhere. Poles on 
the second sheet are called resonances if their real part is 
located above the lowest threshold, and they are called 
virtual states, if they are located below the threshold, but 
on the real axis. It is also possible to have poles on the 
second sheet inside the c omp lex plane with a real part 
lower than the threshold j29|, or on other hidden sheets 
which are often referred to as shadow poles. 

In this study we are interested in branch points on the 
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second sheet in the complex plane, i.e. on the same sheet 
on which the resonance poles are situated. To prove the 
emergence of these branch points, let us start from the 
optical theorem 



= *(27r) 4 £ J d* f 1*(i 



f)T(j^f) (1) 



where T(j — > i) denotes the T-matrix connecting chan- 
nels % and j and d<!> / denotes the phase space of channel /. 
To simplify the argument we assume that the T-matrix is 
in a particular partial wave; below we focus on the singu- 
larities that stem from the unitarity cuts only. Singulari- 
ties like the left-hand cuts, the short nucleon cut [Til l30l| 
or the circular cut, induced by the partial wave projec- 
tion, are ignored in the following for they are irrelevant 
for the argument given. 

To be specific we use the normalization of phase space 
as proposed by the particle data group 31] . Then we 
have for the n-particle phase space 



d$ n (P;pi, ..,Pn) = <5 (4) \ P-^2Pi)l 



d 3 Pi 



(2i:) 3 2E t 



(2) 



where P is the overall center-of-mass (cm.) four- 
momentum. 

To avoid complications, that are irrelevant for the va- 
lidity of the present argumentation, we now focus on 
the diagonal channel i = j. To be concrete we assume 
i = ttN. To further simplify the argument, in addition 
we focus on / = pN as the only relevant intermediate 
ttttN channel. The latter assumption allows us to write 

T(irN -> irirN) = iW(ml„)D(ml„)T(%N pN), (3) 

where D{m\^) denotes the physical p propagator as a 
function of the irir invariant mass and W(m„ n ) is 
the partial wave projected decay vertex, that contains 
also the information on the orbital angular momentum 
I of the decay into tttt. In the following we abbreviate 
to = m^n- 

One can decompose the three-body phase space into 
two subspaces [3l| . 

d® n (P;Pi, -,Pn) = d$j(q;pi, ..,pj) 

x d$ n - j+1 (P; q,p j+ i, ..,p n ){2i:fdm 2 . (4) 



For the example of the p[wTr] N system considered here, 
the first factor d$ refers to the 7T7t phase space in the p 
subsystem at four- momentum q (note that m? = q 2 ), the 
second is the pN phase space at four-momentum P, and 
n = 3, j = 2. With this decomposition, 



d$ j (q;p 1 ,..,p j )\D(m 2 )W(m 2 )\ 2 

= --Im(T>(m 2 )) = p(m 2 ) , (5) 

7T 

where p(m 2 ) denotes the spectral density for the reso- 
nance normalized via 



dm p(m ) = 1 . 



(6) 



We get for the discontinuity of the ttN amplitude from 
the irirN channel 



- (T{itN ->• ttN) - T f (TrN ->• ttN)) 

= (2tt) 7 /" dm 2 p(m 2 ) J d$ 2 (P;g,f> 3 ) 

x \T(ttN ^ pN)(s,m 2 )\ 2 + ... , (7) 

where the ellipses denote contributions from the other 
channels omitted here. The two-body phase space can 
be calculated explicitly. One finds 



d*2(P;q, P3 ) = ^ T dfi , (8) 



with 



p(Vs,to,to 3 ) 
1 



2^~s 



v(s — (^3 + m) 2 )(s — (7713 — to) 2 ) (9) 



for the cm. momentum of the nucleon (particle 3) and 
the pion pair with invariant mass to. 

Using Eq. (J7|) we may thus express the T-matrix 
through a dispersion integral and obtain 



T{ttN -> ttN) 



ds' 



dm p{m ) p(y s' ', m, TO3) / dfl 



\T(ttN -> pN)(s',m 2 )\ 2 



(10) 
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where now the ellipses stand for the unitarity cut con- 
tributions from other channels as well as left-hand cut 
contributions. First of all, there is the three-body cut, 
which drives the inelasticity of the T-matrix. To be con- 
crete, we may write 



p(m 2 ) 



N 
Im 

7T 



1 



r = r 



P 2 t+1 



p 

2t+l 



(11) 



where po is the three-momentum at the nominal res- 
onance mass and N is a normalization factor so that 
Eq. © is fulfilled. The factor (p/po) 2e+1 accounts for 
the centrifugal barrier and p is the pion momentum in 
the p rest frame. Note, p = p{rn,m^,m^) at thresh- 
old, y/s = 27717,- + A/jv, i- e. the p is at rest and the p 
rest frame and overall rest frame coincide. Note also the 
explicit form of Eq. pip is only for illustration. The m- 
dependcncc of the denominator is more complicated in 
general (see, e.g., the Appendix), but the only property 
needed in the following is the presence of poles in the 
spectral function. 

Indeed, the spectral function p(m 2 ) of Eq. con- 
tains a pair of poles located at m 2 = TOq, where too 
denotes the pole position of the p meson, located in the 
complex plane. We may write mo = m p ± iT/2, where T 
denotes the width of the p-meson. 

For the existence of branch points in the complex 
plane, it is sufficient to consider the imaginary part of 
Eq. (|l(jp in the following, or, more correctly, we consider 
the analytic function ST which is ST = Im T for y/s £ R, 
but of course ST / ImT for y/s ^ R (e.g., ST develops an 
imaginary part for complex y/s, whereas ImT does not). 
The function ST can be straightforwardly evaluated, 



ST = - 



dm 2 p(m 2 ) 



4m? 



x p(y/Z,m,M N )p(y/s,m,M N )* L g(y/s,m) (12) 

with p from Eq. @. In Eq. (|12p . we have explicitly de- 
noted a factor of p 2L that comes from the L = 0, 1, • • • 
transition T(irN — > pN). The function g(y/s,m) con- 
tains the integral / dfl over the part of \T\ 2 without these 
centrifugal barrier factors. In general, g{m-\-Mjsi, to) / 0. 
The overall process we consider here as an example is 
shown in Fig. [TJ 

A function f(y/s) has a branch point Zb at y/s = 
Zb, whenever in its integral representation f(y/s) = 

dq f(y/s, q) the function / has a simple pole at q = qo 
and a y/s = Zb exists such that go = a or go = b. For ex- 
ample, the integrand of the two-body phase spac e integra l 
J °° dqq 2 /(y/s - Ei - E 2 + ie), where Ei = y/mf + q 2 , 
has a simple pole at qo = p(s/s,mi,m2) with the on- 
shell momentum p from Eq. Q . Then, the branch point 
is given for the \fs for which qo = (lower integration 
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FIG. 1: The quasi-particle (p) coupling to the stable particle 
N with orbital angular momentum L; the decay of the quasi- 
particle into stable particles (2n) is in i'-wave with respect to 
the quasi-particle cm. frame. 



limit). This is the case for y/s = Zb = m\ + m<x, i.e. the 
branch point is at the two-body threshold. 

With this knowledge, it is straightforward to determine 
the branch points of S T: as discussed before, the simple 
poles of the integrand (spectral function) are located at 
the complex to 2 = TOq which equal the upper integration 
limit of Eq. (fTJ) for y/s = M N + m . 

Thus, without loss of generality, we have shown that 
poles in the spectral function at m = mo lead to branch 
points of the amplitude at the complex scattering energy 
y/s = Mjv + too or 



261,2 



M 



N 



i p ±iT/2 . 



(13) 



More general, the model-independent result is that Zb 
is given by the sum of the mass of the stable particle 
plus to-o, where too is the pole position in the scattering 
amplitude of the subsystem, in this case given by tttt 
which resonates through a p meson. Eq. (fT3"]) has also be 
obtained in Ref . [l7| , starting from an explicit expression 
for the ttttN system, derived from field theory, and in 
which the 7T7t subsystem is boosted. In Appendix [5] we 
will come back to the connection of that formalism to the 
present one. 

The branch points Zb in Eq. (|13|) have been obtained 
by considering the upper integration limit in Eq. (|12[) . 
However, also the lower integration limit can coincide 
with a singularity for a certain y/s: this is the case for 



Zb3 



2to„ 



TO/V 



(14) 



for which the lower integration limit coincides with the 
branch point singularity coming from the factors of p in 
the integrand. The overall analytic structure is shown 
in Fig. [5] The first, physical sheet has the branch point 
Zb3 with an associated cut. If the cut is chosen along the 
real y/s axis like in the figure, the discontinuity of the 
amplitude is given by 2 ST from Eq. (fT2"]) . The branch 
points Zbi and Zbi are in ST, i.e. on the sheet that is 
obtained by analytically continuing the discontinuity of 
the first sheet. They are, thus, on the second sheet, where 
also resonance poles are normally situated. The branch 
points Zb\ and Zb2 induce the new sheets 3 and 4; they are 
analytically connected to the second sheet along the cuts 
induced by Zb\ and Zb2- In Fig. [5] these cuts are chosen 
parallel to the real y/s axis; in Ref. [l7| they are chosen 
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FIG. 2: Analytic structure of the amplitude. There are three 
branch points zti, Z(,2 = ^ti, and zt3- zti and are struc- 
tures in ST and thus on the second sheet. 

parallel to the imaginary y/s axis, which is a convenient 
choice to search for poles. For the numbering of sheets, 
see also Ref. [l7j . 

A. Threshold behavior 
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FIG. 3: Branch point Zt2 in Re ST in the upper \fs half plane, 
for a realistic pN intermediate state. The cut is chosen here 
in the positive Re y^i direction. 



Apart from determining the existence and position of 
branch points, one can also deduce their threshold behav- 
ior, i.e. the functional form of ST close to the three z&. In 
Fig. [TJ the three-body decay is schematically shown. Let 
the quasi-particle (p) couple to the stable particle (N) 
in L-wave in the overall cm. system, while the quasi- 
particle decays into stable particles (2 pions) in ^-wave 
with respect to the quasi-particle cm. frame. 

In the following we will use the explicit form of Eq. I|lip 
to determine the threshold behavior. It is clear, however, 
that the final results do not depend on this particular 
form for the spectral function, but only on the fact that 
the spectral function has poles [right side of Eq. I|15p] and 
the presence of factors of p in Eq. (|12[) that follow from 
the previously given phase space derivation. 

To study the behavior of the amplitude in the complex 
energy plane close to the branch points Zbi,2, complex 
values of m 2 will be needed, and thus the (non-analytic) 
function Im in Eq. (|lip needs to be evaluated to obtain 
a meromorphic expression, 

m — yuio 1 ti\ [rn J 

p(m ) = ^2 : ► — 2 2 • 

n (m 2 — m 2 ) + m 2 T 2 m — m o 

(15) 

The right-hand side shows the behavior of p(m 2 ) close to 
the pole at m = mo; the function h\ does not contain any 
poles or zeros close to mo and thus does not influence the 
threshold behavior. In particular, ft 2e+1 that appears in 
the numerator [cf. Eq. (fTTj) ]. has no zero close to mo and 
can be absorbed in hi. Thus the threshold dependence 
of the branch points Zbi and Zbi does not depend on I , 
which may appear a surprising result. 

To obtain the threshold behavior of the branch points 
Zfci.2 in the complex plane, one inserts the right-hand side 



of Eq. 03) into Eq. |(T2]|. 

ST~ [ dm 2 ^(rn 2 ^ml)^h 2 (m 2 ) 
J m 2 - — TOq 

where we have expanded the argument of the square root 
of the p factors of Eq. (fT2"l) in m 2 , at the point p(-\/s = 
Zb = Mw+niQ, m, M/v) to obtain the power of the leading 
zero from these factors. The function h 2 is again analytic, 
free of zeros close to m = mo, and does not influence the 
threshold behavior. The integral may now be evaluated 
setting this numerator and h 2 constant. The result for 
the threshold behavior of the branch points Zbi,2 is [see 
also Eq. 1(15])] 

ST(z bh2 ) - ~p(v^ J mo,Af Ar ) 2L+1 . 

(17) 

In Fig. [3] we show the branch point Zb 2 in the upper t/s 
half plane [see Fig. [2] for a realistic pN intermediate state 
and L = 0. The branch point is clearly visible, together 
with the cut that in this picture is chosen in the positive 
Re -y/s direction. 

To obtain the threshold behavior for the third branch 
point at z b 3 = 2m n + M/v [see Fig. |2], wc inspect again 
Eq. JT5|). As discussed following Eq. (fTTj) . close to y/s = 
Zb3 = 2m^ + MM [see Eq. (|14|)] the p cm. frame coincides 
with the overall cm. frame, i.e. p = p(m ), and 

thus the ^-wave decay in the p subsystem is also an l- 
wave decay in the overall cm. system. For ^fs in the 
vicinity of Zbs, the denominator of Eq. (|15p is free of zeros; 
however, in contrast to the case of Z6i,2, the numerator 
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f ~ p 2 «+! = p 2 ^ 1 = (m 2 - 4m 2 )(^ +1 )/ 2 does have a 
zero that contributes to the threshold behavior. Inserting 
Eq. (fT5|l (including this factor) in Eq. (Tl2l) and expanding 
the arguments of the square roots of the p factors around 
the zero [cf. Eq. one obtains 



21 + 1 

2t^ /™2 /t™2 \ 2 



(v/i-M«) 

<5T - y dm 2 r (771^- 4771; 
x (m 2 - 4m|) 2 T ti h 3 (m 2 ) , 



(18) 



with a function h 3 free of zeros and poles in the vicinity 
of z b3 . Integration leads now to the threshold behavior 
of z b3 , 



5T(z b3 ) ~ {yfs-{2m„ + M N )) 

~ p(^,A^,2m 7r ) 2£+2L + 4 



(19) 



This corresponds to the opening of the three-body thresh- 
old. Note that even if I = L = 0, the threshold behavior 
is still ~ p 4 , i.e. the standard three-body phase space; 
thus, this threshold opening is always smooth. 



B. The limit of vanishing width 

It is instructive to study the limit of a vanishing width 
of the p-mcson in Eq. ([7]). Then 

p(m 2 ) — > 8{m 2 - m 2 ) for V -> . 
This allows us to perform the m 2 integration to get 

- (T(nN -> ttN) - tUttN -> ttJV)) 

7 

= (2ir) 7 0((^-AI N ) 2 -m 2 p ) J d<S> 2 (P;q,p 3 ) 

2\|2 1 



x |T(7r7V^piV)( S ,m 2 )| 



(20) 



such that Eq. (flUl) reduces to the dispersion integral over 
the standard two-body cut 



T(ttN -> TriV) 



1 



(A/jv+TOp 



= p(Vs',m p ,M N ) 



dfl 



\T(nN^pN)(s',m 2 p )\ 2 



+ 



(21) 



s' — s + te 
The imaginary part which is given by 

ST r ^o = -tt/(VS) p(V5, AfAr) 2L+1 <?(V5, m,) (22) 

has a branch point at yfs — m p + Mjyr, which is simply 
the ordinary two-body threshold on the real y/s axis. As 
r — > 0, the two branch points z/,1.2 in the complex plane 
move towards the real y/s axis until they coincide and 
form this single branch point at yfs = m p + Mjyj. Note 
that there is a factor of V in the numerator of Eq. (|16p . 



but in the limit T — > 0, another factor ~ T appears in the 
denominators from the two poles moving to the real axis, 
that cancels the T of the numerator. Thus, indeed the 
branch point persists in the limit T — > with the result 
given in Eq. 

For the third branch point at z b3 = 2m„+M]yT, Eq. (|18|) 
shows that there are no poles that can prevent the term 
from disappearing in the limit T — > 0; thus, as T — > 0, 
the third branch point fades away. In other words, T — > 
means that the p decouples from mr and thus, in our 
example, the ttttN channel decouples from ttN. 



III. THE RELEVANCE OF BRANCH POINTS 
IN THE COMPLEX PLANE 

As shown in the previous section, whenever there is 
a multi-particle intermediate state with pairwise strong 
correlations, unavoidably branch points show up in the 
complex plane. As we will demonstrate on a particular 
example in this section, their influence on the data might 
well be visible. However, as will be also shown, it is 
in general not possible to deduce the origin of such a 
structure from elastic data only. 

The first model we use is the so-called Jiilich 
model fl4H18l ]. It is a coupled channel meson exchange 
model including the channels irN, r]N, KA, KT, as well 
as 3 effective ttttN channels, namely 7rA, <tN, and pN . 
All these two-pion channels show the mentioned kind of 
branch points [T3|. In Appendix [Al we show the connec- 
tion of the formalism of the Jiilich model to the one of 
the previous section. The Jiilich model allows for a good 
description of the available 7rA^ data in all partial waves 
with j < 3/2 up to an energy of 1.8 GeV and has been 
recently extended to hig her energies, partial waves, and 
additional reactions 



18] 



To be specific we will focus here on the Pll partial 
wave and the region around yfs ~ 1.7 GeV. In this 
energy region, around 300 MeV above the Roper reso- 
nance, signals for another resonance, 7V(1710)Pn, have 



been found in several analyses [31] . It is, however, re- 
markable that in recent analyses of the GWU/SAID 
group 0], there is no sign for this resonance any more. 
Like the GWU/SAID analysis, the Jiilich model contains 
explicitly the branch points z b i^ hi the complex plane at 
y/s = M N + ?77 p ± tX/2 - 1700 ± 75?: MeV. However, 
there are no poles around these energies (the only gen- 
uine pole term in the Pll partial wave is the nucleon, 
while the poles of the Roper resonance are dynamically 
generated [l4j). For the purpose of this study we have 
slightly changed the parameters of the model compared 
to the results of Ref. [l5[ to obtain a good description 
of the GWU/SAID solution. This is shown in Fig. H 
by the dashed lines. The important point here is that 
the theoretical amplitude in the complex plane around 
y/s ~ 1.7 GeV is free of poles, but there is the pN branch 
point. 

To illustrate the difficulties in determining the origin 
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FIG. 4: Fit of the CMB Zagreb model (solid lines) to the 
Pll amplitude provided by the Julich model (dashed lines). 
The "data" points represent the Single Energy Solution of the 
GWU/SAID group @]. 



of structures in the amplitude we fit this Julich model 
amplitude with another model, which does not contain 
the pN branch point in the complex plane. For this, 
we use a Carnegie-Mellon-Berkeley (CMB) type of model 
that has been developed by the Zagreb group @, @, [H, 
[33l | . In this unitary coupled channel model which respects 
analyticity, background plus resonances are provided, but 
all branch points are on the real axis. The result of the fit, 
using two resonance terms, is shown in Fig. [5] by the solid 
lines. As the figure shows, the fit is very precise and, in 
particular, shows no visible discrepancy to the amplitude 
of the Julich model in the energy range shown. 

However, the behavior in the complex plane is quite 
different: as mentioned before, there is no complex 
branch point in the CMB fit by construction; instead, 
a pole is found at 1698 — 130 i MeV which in this case 
might simulate the branch point missing in that model. 

Thus, at a realistic scale of precision, the pN branch 
point does not manifest itself in a unique structure on the 
physical axis; it can be simulated by resonance terms that 
produce poles in the complex plane. Still, the pN branch 
point is a required structure of the 5-matrix, as shown in 
this study, and we have demonstrated that in an analysis 
of partial waves, this and other branch points have to be 
included to avoid false resonance signals, which of course 
can totally distort the spectrum of excited baryonic res- 
onances. 

In such circumstances, one clearly has to consider other 
final states in which the resonance candidate shows a 
clearer signal. As already proposed in Ref. [34|, perform- 
ing global analyses of many different reaction channels 
within one theoretical ansatz is a much cleaner way to 
determine the resonance spectrum than increasing the 
precision of a partial wave for one reaction. 

First steps within the coupled-channel Julich model 
have been undertaken in this direction through the in- 



clusion of some pN data [l4|, rjN data [T^], and, most 
recently, K + T, + data [ll]- For the isospin I = 1/2 sector, 
we expect the inclusion of K°A data to further clarify the 
role of the N(17W)P n [see also Ref. |U]. 

Thus, the aim of the present short exercise is not to 
discard the existence of the much-debated AT(1710)Pn 
as such; rather, we have shown that branch points in the 
complex plane are relevant; in their absence, resonances 
may be needed to simulate them, and, thus, the extracted 
baryon spectrum can be easily distorted. 



IV. CONCLUSIONS 

Using only general properties of the S'-matrix we have 
shown the existence and determined the position of three 
branch points induced by intermediate quasi-two body 
states. Those are three-body states in which two particles 
are so strongly correlated that the scattering amplitude 
of this subsystem has a pole. A pole in the subsystem 
necessarily leads to the appearance of branch points in 
the complex y/s plane of the overall ttN amplitude. This 
result is model-independent because it does not depend 
on any particular parameterization, but only on analytic- 
ity and general properties of the three-body phase space. 
We have also determined the threshold behavior of all 
branch points, which depends on the orbital angular mo- 
menta of the two decay processes involved. Finally, on 
the example of the Pll partial wave, it has been shown 
that branch points in the complex plane are relevant in 
partial wave analysis: if a theoretical amplitude does not 
contain the branch points, false resonance signals may 
be obtained. To allow for a reliable extraction of the 
baryon spectrum, it is thus mandatory to include also 
these branch points in the analysis. 
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Appendix A: Spectral representation of the Julich 
model 

In this Appendix the connection of the field theoretical 
formalism, used in the Julich model of hadron exchange, 
to the formalism used in this study is outlined, up to 
overall normalization factors. For further details of the 
formalism used in the Julich model, we refer to Ref. [l7| • 
For the example of the pN propagator, that is considered 
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here, the propagator on the real axis is given by 



g P N(Vs,k) 



v/5 - E N (k) - E°Jk) - S(z p (Vi, k), k) 



(Al) 

where En is the nucleon energy, E° is the p energy using 
the bare p mass and E is the p self energy, where z p (y/s, k) 
is the boosted energy for the p subsystem. The explicit 
form of z p (v / s, k) is quoted in Ref. [17} but for the present 
discussion the only needed property is that z p (yfs,k = 
0) = ^fs — AIn- The propagator g p N is iterated in the 
multichannel scattering equation, but to investigate the 
analytic structure it is sufficient to consider the one-loop 
amplitude 



G pN {^~s) = j dkk 2 g pN {^~s,k) (A2) 
o 

where for simplicity we have omitted the form factors 
that regularize this divergent expression. One can rewrite 
the Dyson-Schwinger representation of Eq. (|Alj) with the 
spectral function 



S(uj, k) = Im ■ 



- E N {k) - E^{k) - Y,(z p (uj, k), k) 



resulting in the Lehmann representation 

oo 

S{oj,k) 



9 p n{Vs, k) 



duj ■ 



fs — 10 + ie 



(A3) 



(A4) 



2m.T,+M N 



For the imaginary part of the pN loop G p n{\/s), one 
obtains: 



Im G p n{Vs) = Im J dkk 2 g pN (y/s,k) 
o 

oo k± 

= -n J dkk 2 S{^~s,k) = -7T J dkk 2 S(^,k) . (A5) 



The last equality shows that the integration can be cut 
at k = k\ as for k > k\ the spectral function is zero be- 
cause then z p (^/s 1 k) < 2m n . In particular, k\ is given 
by z p (y/~s, k\) = 2m- !T . Note that the explicit evaluation 
of the integration limits as done here is necessary if one 
wants to use the spectral representation in the complex 
yfs plane. This has been shown recently in the context of 
Feynman parameterized loops [3(| : the integration limits 
have to be analytically continued for complex ^fs to ob- 
tain the analytic continuation of the loop itself, and for 
this they need to be known explicitly 



Eq. (|A5|) can be rewritten as 



Im G pN (Vs) = 



dm 



S(^,k on (m))m 



X (-7T) 



k on {m)E° n E° 



(A6) 



with 



fc on (m) = p(^,m,m„), E° n = V 'ml + (fc°») 2 , 

E% = \/™ 2 + (^) 2 . (A7) 

and p from Eq. @. The lower integration limit mi is 
given as the solution of z p (y / s, fc on (mi)) = 2m 7r . The 
second fraction in Eq. (|A6|) can be compared to the imag- 
inary part of the well-known [TtJ propagator of two stable 
particles M and N, 



Im G 



stable 



k° n (m = y/s)EffiEff 



(A8) 



Thus, the imaginary part of a loop with one stable and 
one unstable particle can be expressed as an integral over 
a distribution of imaginary parts of the form of Eq. (|A8|) . 
Comparing Eq. (|A6|) to Eq. p2[) . one sees the formal 
similarity: there is an integral of a spectral function, 
that has poles [cf. Eq. (|A3|) ]. together with the factor 
k on (m) = p(y/s, m^m^), and both ingredients produce 
the three branch points Zf)i,2,3 as has been shown in the 
main text (we have omitted here the additional 2L pow- 
ers of p for simplicity) . There is a difference in the cho- 
sen parameterization in terms of the spectral function 
[compare p(y / s, m, m^) in Eq. (|A7[) vs. p{^fs,m,Mjq) in 
Eq. (|12p], but this does not change the position of the 
branch points. 

Indeed, k on = for the upper integration limit m ~ 
^fs — m-K and thus z p (y / s, 0) = ^fs — Alfq. The poles of the 
p resonance in the spectral function S are located at the 



complex z. 



and consequently the integration limit 



equals the pole position for ^fs = Z(,i,2 = M 



N ■ 



which 

is indeed Eq. fT3J). The singularity at k° n = 0, coming 
from the factor k on (m) in Eq. (|A6|) . is also reached if 
to = t/s — m n . It is easy to show that this to equals the 
lower integration limit toi for y/s = 2m 1T + M/v and thus 
Eq. (|A6[) indeed provides also the third branch point Zb3 
from Eq. QMI) . 

In Ref. [17I ]. the amplitude has been analytically con- 
tinued to the complex plane using contour deformation. 
In fact, one could use the representation of Eq. (|A6[) for 
the same purpose in principle. As shown in the main text, 
for this, one has to respect the analyticity of the spectral 
function, i.e. the Im () function in Eq. (|A3[) needs to be 
explicitly evaluated like in Eq. (f*H))) . Second, and this is 
an additional complication, the self energy £ itself has a 
two-sheet structure and the corresponding cut needs to 
be rotated as specified in Ref. [TtJ • This cut in £ induces 
the cut of branch point z^ in the overall ttN amplitude. 
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Apart from this and a carefully chosen integration path complications, and the spectral representation allows for 
for the m- integration of Eq. (|A6[) . there are no additional an alternative way of analytic continuation. 
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